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Abstract. Various extensions of the Standard Model predict the existence of hidden 
photons kinetically mixing with the ordinary photon. This mixing leads to oscillations 
between photons and hidden photons, analogous to the observed oscillations between 
different neutrino flavors. In this context, we derive new bounds on the photon-hidden 
photon mixing parameters using the high precision cosmic microwave background spec- 
tral data collected by the Far Infrared Absolute Spectrophotometer instrument on 
board of the Cosmic Background Explorer. Requiring the distortions of the CMB 
induced by the photon-hidden photon mixing to be smaller than experimental upper 
limits, this leads to a bound on the mixing angle xo ^ 10^^ — 10^^ for hidden photon 
masses between 10~^'*eV and 10^^ eV. This low-mass and low-mixing region of the 
hidden photon parameter space was previously unconstrained. 
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1. Introduction 

The last decade has witnessed the blooming of modern cosmology, supported by a 
consolidated model of particle physics and a huge amount of new observational data 
with unprecedented precision (see for example the latest WMAP team results [1,2]). 
Every observation has been shown to fit with the predictions of the so-called A-CDM 
model, a spatially flat Friedman universe, with ~ 10~® baryons per photon and whose 
energy composition is now dominated by two unknown ingredients: dark matter and 
dark energy. At this stage of concordance, observational cosmology can be used to test 
the existence of non-standard particle physics. 

One of the most powerful probes used to constrain exotic physics is represented by 
the observations of the cosmic microwave background radiation. Standard cosmology 
predicts relic radiation from the big bang which, being originally in thermal equilibrium 
with matter, suffered a soft decoupling and, therefore, still features today a perfect black 
body spectrum. As it turns out, during most of the history of the universe this perfect 
blackbody has been completely unprotected against distortions. Therefore, any new 
particle beyond the standard model which can interact with photons can potentially 
distort the blackbody spectrum, providing us with a glimpse of the existence of such a 
particle. 

These distortions are expected to be more severe if the new particles have feeble 
couplings, so that their interactions with photons do not allow to establish thermal 
equilibrium (but of course strong enough to produce some effect), and if they have 
small masses, such that photons can annihilate into them and/or scattering processes 
are not suppressed by a heavy mass scale. Since the cosmic radiation spectrum freezes 
out when the temperature of the universe falls below the ~ eV scale, the spectrum of the 
cosmic microwave background (CMB) can, therefore, be an excellent probe for weakly 
interacting sub-eV particles (WISPs). 

Among the most elusive WISPs we find low mass hidden photons (HPs), gauge 
bosons of a W(l)jjjj gauge symmetry having kinetic mixing with the ordinary photon. 
Due to the mixing one expects photon-HP oscillations driven by the mass difference of 
the two particles. Clearly, as this mass becomes smaller the oscillation length grows 
and oscillations are harder to detect. In this respect, the CMB provides an excellent 
probe for photon oscillations into low mass HPs, since the beam-line is the longest at 
our disposal, namely the whole universe. Indeed, hidden photons were the first WISP 
candidates to be confronted with CMB data in the early 80 's [3]. 

Recently, new intriguing ideas and experimental techniques have been proposed 
to achieve possible detections of these still elusive particles [4-6]. Moreover, the 
astrophysical and cosmological role of HPs is currently under exploration [7-10]. 

The purpose of this paper is to revisit the CMB bounds on photon-HP mixing using 
the most precise available observations of the CMB spectrum, those provided by the 
Far Infrared Absolute Spectrophotometer (FIRAS) on board of the Cosmic Background 
Explorer (COBE) [11, 12]. The high precision of this measurement, which confirmed 
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the blackbody nature of the spectrum at better than 1 part in 10^, has been aheady 
exploited to constrain hidden photons [13] and also other WlSPs such as axions [14], 
or radiative neutrino decays [15] or millicharged particles [16]. However, the already 
existing constraints [3,13] did not properly take into account the refractive properties 
of the primordial plasma and are thus incomplete. The medium effects are especially 
important when we realize that a resonant photon-HP conversion is possible, similar 
to the Mikheev-Smirnov-Wolfenstein (MSW) effect in the neutrino case [17-19]. This 
resonant conversion is much stronger than the vacuum oscillations considered so far and, 
therefore, the bounds we obtain are much stronger. 

The plan of our work is as follows. In Section 2 we summarize the relevant formalism 
concerning the photon-HP mixing. In particular, we describe how the medium effects 
modify the oscillations in our system. In Section 3 we present our analytical recipe 
to calculate the photon-HP conversion probability across a resonance in the expanding 
universe. In Section 4 we describe our simplified model for the effective photon mass, 
induced by the primordial plasma. In Section 5 we describe the constraints coming 
from spectral CMB distortions. In Section 6 we outline in a general way the strength of 
photon-HP mixing bounds that can be obtained from an observed modification of the 
photon flux from a generic astrophysical source. Finally, in Section 7 we compare our 
new cosmological bound with the other ones existing in the literature and we draw our 
conclusions. 



2. Photon-Hidden Photon Mixing 

A hidden photon is the gauge boson of a gauge ^/(l)hid symmetry under which all the 
SM fields are uncharged, and thus remains hidden in our world. At energies above the 
electroweak scale, particles charged under W(l)jjj^ and the SM hypercharge are likely 
to exist and act as messengers between the two "sectors" of low-energy physics. Their 
effects will produce effective operators in the low-energy theory that contain the SM 
and the hidden photon. Typically these operators are irrelevant, suppressed by the 
heavy particle masses. But there is still one marginal operator, the so called kinetic 
mixing [20] whose effects are not in principle suppressed by heavy masses and thus can 
lead to relatively large effects. The low energy Lagrangian of such an extension of the 
SM would be then [8] 

£ = - - F F^" --B B^" 

+ '^B,„F'- + '^n.l.B.B^ + jS^A, , (1) 

where A^, 5^ are, respectively, the photon and the hidden photon fields, and B^i, 
their respective field strengths and j'^^^ is the electromagnetic current. We have included 
a hidden photon mass my which we will treat as a free parameter. Typical predicted 
values for the mixing angle Xo in realistic string compactifications range between 10~^^ 
and 10-2 [21] (see also [22,23]). 
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The kinetic mixing term can be removed, leading to a canonical form of the kinetic 
lagrangian, by the following change of basis 

{A,B}^{An,S} , (2) 

where 

Ar = cos xoA , (3) 

S =B-smxoA . (4) 

Let us call 7,7^ the quanta of the Ar and S fields, respectively. The basis {7, 7s} 
can be called the "interaction basis", since the photon 7 is the state that interacts 
with SM charged particles and 7^, being orthogonal to 7, is completely sterile. In the 
interaction basis the kinetic lagrangian is diagonal but the mixing angle Xo appears in 
an off-diagonal term in the mass-squared matrix, 

^2^/ m2,sinXo "^y sin xo cos Xo \ 
\ m^, sin Xo cos Xo rn^' cos^ Xo / 

The mass matrix can be diagonalized through the unitary matrix U 

^^/cosxo -sinxo\ 
V smxo cosxo / 
that allows to identify the two "propagation states" 

where 71 is mostly photon-like and massless, while 72 has mass my and is close to the 
sterile state. The mismatch between the interaction {7,7^} and propagation {71,72} 
states by the mixing angle Xo is well-known to produce 7 — > 7s oscillations [24], with a 
conversion probability in vacuum given by 

P^^^, = sin2 2xosin^(myL/4tu) , (8) 
where L is the path length and ut is the photon energy. 

In analogy to the neutrino case [25], the 7 7s oscillations are modified by the 
refraction properties of the medium. In the primordial plasma, photons acquire a non- 
trivial dispersion relation which can be parametrized by adding an effective photon 
mass to the Lagrangian. This mass depends on the properties of the medium, as 
we will show explicitly in Section 4. This is generally complex, reflecting the absorption 
properties of the plasma, but in the context relevant for the present paper the imaginary 
part is negligible. In this case, the effective mixing angle is related to the vacuum by [7, 8] 

sin 2xo 




sin2x 



sin^ 2xo + (cos2xo-0' 



COS2X = , . "^''■•'-^ ^ . (9) 



sin 2xo + (cos2xo - 



where the parameter 



e = m^/m^, (10) 
measures the significance of the medium effects. 
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3. Photon-Hidden Photon Oscillations in the Expanding Universe 

To obtain the bound on photon-hidden photon mixing, we have to compute the fraction 
of CMB photons that would oscillate into invisible 7^. The photon effective mass squared 

is generally proportional to the density of charged particles in the medium, so it 
relaxes as the photons propagate in the expanding primordial plasma. We are then 
facing a complicated problem of oscillations in an inhomogeneous medium. 

For sufficiently early times, ^ niy and, therefore, the photons are very close 
to be both interaction and propagation states at a given time. In this case, oscillations 
into 7s are suppressed, 

,^^1, x~^^/2 (medium suppression) . (11) 

As the universe expands, ^ eventually reaches 

C, = cosxo ) X ""^ (resonance condition) . (12) 

When this condition is fulfilled, resonant photon-hidden photon conversions are possible, 
analogous to the well-known Mikheev-Smirnov-Wolfenstein (MSW) effect in the neutrino 
case [17-19]. As ,^ decreases below unity, 

^ <^ 1 , X ^ Xo (vacuum oscillations) (13) 

7 7s oscillations take place as if they occurred in vacuum. 

In vacuum, oscillations are proportional to sin^ 2xo [see Eq. (8)] and since the 
blackbody nature of the CMB is established experimentally to an accuracy of ~ 10~^ 
we can easily exclude mixing angles ^10~^. However, during a resonance sin 2% ~ 1 
and the conversion 7 — * 7s can be much stronger. In this paper we want to focus on 
the range of HP masses that can undergo such a resonant transition, and are, therefore, 
strongly constrained. 

When the photon production and detection points are separated by many oscillation 
lengths, then on both sides of a resonance, = m^/, the oscillation patterns wash out 
and the transition probability is given by [26] 

P-/->-,s = (p^^)^^^ cos 2x , (14) 

where xo, X ^i-re the mixing angles at the detection and production points considered to 
be in vacuum and high density, respectively, and p is the level crossing probability. This 
latter takes into account the deviation from adiabaticity of photon-HP oscillations in 
the resonance region. In particular, one has p = for a completely adiabatic transition 
and p = 1 for an extremely nonadiabatic one. 

As we will see, we are going to bound mixing angles much smaller than xo ~ 10~^, 
so for simplicity we can already take cos 2xo — — cos 2% = 1 in Eq. (14) which, therefore, 
takes the much simpler form 

P,^,^c^l-p. (15) 

The FIRAS sensitivity will allow us to bound 1—p < 10~^, requiring thus a non-adiabatic 
resonance. 
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The crossing probability p for photon-hidden photon resonant conversions can be 
obtained using the Landau-Zener expression 

p ~ exp(— 27rrA; sin^ Xo) ? (16) 
where k = m'^,/2u} is the 7 — > 7s vacuum oscillation wavenumber and 

(17) 

is a scale parameter to be evaluated at the location where a resonance occurs. 

This expression has been widely used in solar and supernova neutrino 
oscillations [25]. It represents an accurate ansatz to calculate the level crossing 
probability for general density profiles (see, e.g. [27]). In particular, we observe that it 
reproduces the correct limits p ^ for r ^ 00 (adiabatic limit) and p ^ 1 for r — * 
(extreme non-adiabatic limit). 

Finally, under the approximations outlined before, we can approximate the non- 
adiabatic conversion probability as 

-1 

(18) 

t — £rcs 

The half- width width of the resonance is, according to Eq. (9), 5^{t) ~ sin2xo; which 
corresponds to a time scale 

Tr ~ r sin 2xo • (19) 

Due to the smallness of the vacuum mixing angle, the resonance is very narrow, so that 
it is reasonable to take into account the deviation from the adiabaticity only at the 
crossing point. 



rfln mUt) 



dt 



2'Krkx, 



UJ 



d\n mt^{t) 
Jt 



4. Cosmological profile 

From Eq. (18) it is clear that the only information needed to calculate the conversion 
probability is the profile of the photon effective mass along the cosmological line of 
sight. The primordial plasma is thought to be electrically neutral and composed mainly 
of hydrogen and helium in a fraction per mass Yp = m-nc/fnn ^ 0.25. The effective mass 
squared has a positive and a negative contribution [28] from scattering off free electrons 
and off neutral atoms J 

^ - 2uj^{n - 1)h 

^ 1.4 X 10-2^ (x, - 7.3 X 10-^ (^)' (1 - X,)^ £^ eV^ , (20) 

where = ATrane/mg is the plasma frequency with a the fine structure constant, the 
electron mass and rig the free electron density. We have written the ionized fraction of 
hydrogen as = Ue/up with Up the proton density. The indices of refraction of neutral 

X The vacuum magnetic birefringence [29] due to a primordial magnetic field is negligible [30]. 
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hydrogen and helium are (n — 1)h = 13.6 x 10~^ and {n — l)uc = 3.48 x 10~^ in normal 
conditions [28], rather insensitive to u. Since the fraction of electrons corresponding to 
helium is ~ 13% and (n — 1)hc {n — 1)h we have neglected the effects of helium in 
Eq. (20). 

The photon frequency u and the proton density Up are given in terms of redshift z, 
the photon energy today c^o, the CMB temperature today Tq and the baryon to photon 
ratio r/ ~ 6.7 X 10"^° as [31] 

CO = uoil + ^) ; = (l - ^) ^^^0(1 + ■ (21) 

Note that for a negative effective mass squared a resonance is not possible. However, 
the negative contribution is proportional to u;^(l — X^) and thus is unimportant for 
sufficiently small frequencies (small or late times when the redshift z is small) and/or 
large ionization fractions, ~ 1. Interestingly, the smallest frequencies are the most 
precisely determined by FIRAS. 

Since the conversion probability Eq. (18) is proportional to r and u!~^ which both 
grow with decreasing redshift, for a fixed HP mass the later the resonance the more 
adiabatic the transition will be and the stronger our constraints on the mixing parameter 
Xo will be. 

The history of the ionization fraction Xe{z) is extremely complex. Above a 
temperature T ~ 0.5 eV (redshift z ~ 1100) hydrogen is fully ionized. As the universe 
temperature decreases, photons cannot ionize hydrogen efficiently and electrons and 
protons slowly combine. This makes the universe very transparent to radiation, indeed 
releasing the photon bath which we see today as the CMB. This epoch of so-called 
recombination has been studied in great detail in [32] from which we can take the values 
of Xe as a function of redshift, shown in Fig. 1. Later on, the universe becomes ionized 
again due to ultraviolet radiation from the first quasars or population III stars. The 5 
year data of the WMAP mission constrains the redshift of an instantaneous reionization 
to 2; = 11.0 ± 1.4 with 68% confidence level by CMB polarization studies [1,2]. At 
the same time reionization can be studied with the spectra of high redshift quasars, 
which show so-called Gunn-Peterson troughs due to absorption of light at Lyman-a 
frequencies. The study of several quasars by the Sloan Digital Sky Survey showed that 
reionization should end around a redshift 2 ~ 6, and, therefore, should be an extended 
process taking place between z ~ 6 — 11 [33]. 

In Fig. 1 we show a possible profile for the cosmological history of X^ and TTly clS db 
function of redshift z. We can learn several important things from this figure. Let us 
focus on the right panel. We can first compare the thin line, which corresponds to a 
cosmology when atoms are always ionized, i.e. X^ = 1, with the more realistic colored 
lines which includes the ionization history of the left panel and assumes uj/T = 1, 3, 4, 10, 
respectively. We see that during the dark ages (6 < z < 1000), Xg drops very much and 
the resonances are moved to higher redshifts compared with the X^ = 1 case. As we 
commented before, the later a resonance the stronger it is, and, therefore, the inclusion 
of the ionization history tends to decrease the bounds on xo- 
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z 



Figure 1. Left panel: The ionization fraction of hydrogen. For redshifts above z = 11 
the hne is taken from the data of [32]. The re-ionization period between z = 6 and 
z = 11 is modeled. Right panel: Effective photon mass as function of redshift. The 
thin line uses X,, ~ 1, the blue, magenta, green and red lines are for lu/T = 1, 3, 4, 10, 
respectively. The two sharp dips bound the region where m-y becomes imaginary. 



On the one hand, due to the smallness of Xg during the dark ages, the index of 
refraction of neutral hydrogen can dominate over the contribution from free electrons 
turning the effective photon mass squared negative (see Eq. (20)). Since this negative 
contribution is proportional to the photon frequency squared, the region of redshifts 
for which this happens is broader for higher frequencies. From Fig. 1 we see that this 
does not happen for u;/T^3 but for cj/T = 4 it forbids a resonance in the range of 
redshifts 400^^:^700 and for uj/T = 10 the range without resonance is enlarged to 
100<z<1000. 

Let us recall that this does not forbid a resonant transition but makes it happen 
at higher redshifts than those which correspond to u/T ~ 1. Consider as an example 
the case my = 4 x 10~^^ eV. For frequencies uj/T <^2 the resonance, = my, occurs 
at z ^ 400, for uj/T = 4 at z ~ 700 and for cu/T = 10 at z ~ 1000. Moreover, for 
these delayed resonances it becomes clear that the value of r, which is related to the 
inverse of the derivative of m^ is smaller than for resonances for low u/T. Because 
of these two facts, the resonances for high u/T photons during the dark ages will be 
much less adiabatic than those for lower u/T and, therefore, they will play a minor 
role in our bounds. There is only one possible exception to this conclusion. For very 
small masses my^lO"^^ eV there is no resonance at all for low u/T but still there 
are resonances for uj/T^A during the dark ages. We have checked however that such 
resonances are extremely non-adiabatic and no interesting bounds can be derived from 
them. Therefore, in the following discussion we will not discuss resonances for u/T^4, 
although we will include them in our bounds. 

Regarding the resonant conversion 7 — 7s, and focusing on uj/T<,3, we see three 
different possibilities in the homogeneous universe: 
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• fT^Y ^2 X 10 eV: there is only one resonance crossing 

• 2 X 10"^^ eV ^my ^10^^^ eV: there are three level crossings respectively after, 
during and before reionization. 

• my <, 10~^^ eV: no resonance is possible. 

Our analytical ansatz for the crossing probability [Eq. (16)] has been shown to 
reproduce accurately the transitions also for nonmonotonic density profiles, where 
multiple crossings could arise (see, e.g., the supernova neutrino case [27]). In this case, 
assuming the factorization of n different crossings, one obtains [27] 



where the conversion probability for a resonance is still given by Eq. (18). As expected, 
since we know that photon depletion has to be very small, (and, therefore, also a possible 
re-creation) the dissappearence probabilities at different crossings just add up. 

As commented previously, in the presence of multiple crossings (2 x 10~^^ eV 
^my^lO"^^ eV), the relevant crossing is the one occuring after reionization (at 
z < 6), since the other two will be less adiabatic. This situation further simplifies 
our calculations. Since the crossings during reionization do not play a role, we need 
to evaluate them only for z <6 (for my < 2 x 10"^^ eV) and for, say, z>70 (for 
my > 2 X 10^13 eV). 

During the epoch of matter domination, the small density inhomogeneities already 
present at CMB decoupling grow and the density profile along a line of sight can be 
quite complicated at the smallest redshifts, with lots of resonances at galactic length 
scales. However, according to Eq. (23) further resonances tend to increase the transition 
probability, thus tending to strengthen the constraints. It is, therefore, conservative to 
restrict oneself to the most adiabatic resonances in the smooth background, represented 
in Fig. 1. 

Furthermore, the optical depth r for crossing an object of typical size / with 
number density n along a cosmological line of sight of length s ~ 4Gpc is r ~ nPs ~ 
0.004(10~^/Mpc~'^)(//10 kpc)^ which is small for galactic objects. Within an order of 
magnitude this is consistent with galaxy number counts which yield ~ 10® sr~^ [34]: 
Given that one galaxy subtends a solid angle Q ~ (10 kpc/10'^ Mpc)^ sr ~ 10"^" sr, 
this gives r ~ 0.01. Similarly, counts of galaxy cluster with mass >10^^Mq yield 
~ 10~^Mpc~'^ [35], or ~ 10^ sr"^. With a typical solid angle subtended of f2 ~ 
(1 Mpc/lO"^ Mpc)^ sr ~ 10"^ sr, this also gives r ~ 0.01 for the average number of galaxy 
clusters crossed by a given line of sight. 



p=o 1- n (i-2p.) 



(22) 




we can write 



(23) 
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4.I. Adiabaticity parameter 

Given our density profile it is now straightforward to compute tlie adiabaticity parameter 
r for each crossing and from it, the photon disappearance probabihty. For masses 
my < 10"^^ eV this would be very much dependent on the inhomogeneities at small 
redshifts, and we prefer not to treat this case. We are left with HP masses above the 
average effective cosmo logical plasma mass of photons at zero redshift. 

Recall that in general r oc t so the most adiabatic crossing is the latest, when the 
universe expands slower and, therefore, sweeps smaller ranges of in a given time. 
This argument also allows us to discard the effect of spatial inhomogeneities unless they 
have sizes comparable with the size of the universe. It also makes us neglect a possible 
back-reaction during the reionization crossing. 

Let us then write, 

d log ml^ _d log dz 

Jt ~ Tz ' ^ ' 

with 

dz I 

— = -ifo(l + 2) V^A + ^m{l + Zf + a.(l + ZY , (25) 

where we fix Hq = 70 km s~^ Mpc~^, Qm = 0.27, Qx = 0.73, fi^ = 6 x 10~^ consistent 
with the recent determinations of cosmological parameters [33]. 

The function dX^/dz can be computed numerically. However, since the crossing 
during reionization does not play a role, we need this function only for z <^Q (for 
my < 2 X 10~^^ eV) and for 2;>70 (for my > 2 x 10~^^ eV). Below 2; = 6 we can 
take it to zero, and for 2 > 70 we can obtain it from the following fitting function 

log^„Xe = -3.15 ^ , (26) 
e 160 -|- 1 

valid at the % level. 
5. FIRAS Bounds 

The CMB spectrum measured by FIRAS fits extremely well to a black body spectrum 
at a temperature To = 2.725 ± 0.002 [12]. The energy range of the CMB spectrum 
measured by FIRAS [11] is 2.84 x 10~'^eV < ooq < 2.65 x 10"^ eV, corresponding 
to 1.2 <a;o/^o ^ 11-3. In that region, the CMB blackbody becomes unprotected to 
distortions below a cosmic temperature ~ keV, which corresponds to a photon mass of 
~ 10""^ eV. On the other hand, today the average plasma mass for photons is as low as 
2 X 10~^^ eV. If HPs exist with a mass between these two values they will be produced 
resonantly and leave their imprint on the CMB. 

The accuracy of FIRAS constraints P-y-,^^ ^ / — 10~^, which, using Eq. (18) leads 
to the bound 
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Figure 2. Bounds from distortions of the CMB blackbody due to 7 ^ 7s photon 
depletion: 95% C.L. (sohd) and 99% C.L. (dashed). 



In order to sharpen this bound, we have considered the distortion of the overall 
blackbody spectrum. 

To this end we use the COBE-FIRAS data for the experimentally measured 
spectrum, corrected for foregrounds [11]. Note that the new cahbration of FIRAS [12] 
is within the old errors and would not change any of our conclusions. The = 43 data 
points <l>°^'' at different frequencies Ui are obtained by summing the best-fit blackbody 
spectrum to the residuals reported in Ref. [11]. The errors a^^^ are also available. In the 
presence of photon-hidden photon conversion, the original intensity of the "theoretical 
blackbody" at temperature T 

$°(^,T) = ^[exp(.;/T)-l]-^ (28) 

would be deformed to ^{uj, T, xo, ^7') = T)[l — P^^^^ui, Xo? m^i)]. We then build 

the reduced chi-squared function 



cxp 



(29) 



We minimize this function with respect to T for each point in the parameter space 
A = (m^/,xo)5 T is an empirical parameter determined by the minimization for 
each A rather than being fixed at the standard value Tq = 2.725 ± 0.002 K. 

In Fig. 2 we show our exclusion contour in the plane of m^i and xo- The region 
above the continuous curve is the excluded region at 95% C.L., i.e. in this region the 
chance probability to experimentally obtain larger values of is lower than 5%. We 
also show the corresponding 99% C.L. contour which is very close to the 95% contour so 
that another regression method and/or exclusion criterion would not change the results 
very much. Note that for masses m^i for which the resonant conversion takes place 
during the dark ages such a resonance is only possible for lu/T < 2 for which the positive 
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contribution to dominates over the negative one. Only the lowest four points in the 
FIRAS data set satisfy u/T < 2 and are thus distorted. On the other hand, for larger 
(my > 10~^° eV, cf. Fig. 1) or smaller masses (2 x 10~^^ eV < my < 2 x 10~^^ eV), 
the resonance can happen up to Ci;/T = 10 and we can use the whole FIRAS data set. 
This explains why the bound we obtain is stronger in this latter region. 

Finally, we comment that in our analysis we have assumed that hidden photons are 
produced only by oscillations. In this sense, we have neglected a primordial population of 
hidden photons, that could have subsequently been converted into photons, producing 
an additional distortion of the CMB. Our assumption is reasonable if we neglect the 
presence of primordial charged hidden particles, by whose annihilations hidden photons 
could have been produced. 

6. Astrophysical bounds 

The constraint of Eq. (27) holds for general astrophysical sources whose photon flux at 
frequency u is known to deviate less than a fraction / from a model prediction in the 
absence of mixing with hidden photons. 

As long as the contribution of neutrals to m^ is negligible, the resonance condition 
is m^ = uj"^ . The resonance thus occurs at the plasma density 



Such resonances can directly influence the photon flux only if the optical depth at 
the resonance is smaller than unity. Since the optical depth ^ c-pne r with — 
6.7 X 10~^^ cm^ the Thomson cross section, from Eq. (30) we obtain the condition 



for the scale over which the mixing potential varies. Since furthermore only photons 
with frequency u'^Up can propagate, we conclude from Eq. (27) that independent of 
the details of the astrophysical system, only mixing parameters satisfying 



can be constrained by observing photon fluxes from astrophysical or cosmological 
objects. Note for example that for niy ~ meV, the best possible bound from 
direct observations is thus x 10"^/^^^. Any mixing parameters smaller than 

this, therefore, has to be constrained by indirect methods (e.g. stellar cooling) or 
experimentally. This theoretical best possible bound from astrophysical sources is 
plotted for / = 1 in Fig. 3. Our CMB based bounds are weaker than this best possible 
bound essentially because the scale r in Eq. (27) over which the cosmological plasma 
mass varies is considerably smaller than the photon mean free path Eq. (31) that would 
lead to the most adiabatic resonance possible for unabsorbed photons. 




(30) 




(31) 




(32) 
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Figure 3. Bounds from resonant 7 ^ 7s depletion of tlie CMB blaclcbody as 
constrained by FIRAS data in this work (gray region) and in [8] (black region). The 
gray diagonal line separates the region where the resonance happens at small damping 
(left, this work) or at strong damping (right, cf. [8]). Plotted for comparison are bounds 
from tests of the Coulomb l/r^ law [36,37], magnetic fields of Jupiter and earth [38], 
photon-regeneration-experiments [39-44] , arguments of the lifetime of the Sun and the 
CAST search of solar axions [7,45]. The solid black line indicates the best possible 
bound Eq. (32) that can be obtained from astrophysical or cosmological sources whose 
photon flux is known to be unmodified by photon-HP mixing to order unity. 

7. Conclusions 

In this paper, we have revisited the bounds on photon-hidden photon oscillations coming 
from cosmology, deriving updated constraints from the high precision CMB spectrum 
data collected by the FIRAS instrument on board of COBE. We have also commented 
about complementary bounds from a depletion of the photon flux from astrophysical 
sources. Previous studies [3,13] were derived in the pre-COBE era and/or more 
importantly they lacked a detailed treatment of the effects of the plasma medium on the 
photon-hidden photon oscillations. This has motivated us to re-evaluate the bounds. 
This problem has presented interesting analogies to similar situations encountered in 
neutrino oscillation physics. In this regard, we had the benefit to apply to our case 
different analytical recipes developed in neutrino oscillation studies. 



Microwave Background Constraints on Mixing of Photons with Hidden Photons 14 

The result of our analysis leads to a bound on the mixing angle Xo ^ 10^'' — 10~^ 
for hidden photon masses between 10~^^eV and 10~^eV. In Fig. 3 we plot our FIRAS 
bound together with the other ones existing in the literature (see [46] for a complete 
review). It turns out that our new cosmological bound excludes a region of low-mass 
and low-mixing angle in the hidden photon parameter space, that was unconstrained 
by previous arguments. As a result of our new bound, it is unlikely that hidden 
photons with masses smaller than 10~^eV can play a cosmological role. Conversely, for 
meV masses, resonant photon-hidden photon oscillations happen after nucleosynthesis 
but before CMB decoupling, increasing the effective number of neutrinos but also the 
baryon to photon ratio with interesting cosmological consequences [8]. The mixing 
angles required for this effect could be probed in current laboratory experiments. 

We comment that our cosmological bound on hidden photons is based on an 
approximate treatment of the plasma environment. This would leave open possible 
improvements of our limit when a more accurate description of the primordial plasma 
will be achieved. Moreover, further improvements of our cosmological bound for 
extremely small masses and mixing could be reached observing a possible depletion 
of photons from astrophysical sources, as we proposed in Sec. 6. In this context, for 
example it still remains to investigate if photon-hidden photon oscillations could have 
some impact on the apparent dimming of supernovae la, mimicking the acceleration of 
the universe. 
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